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ABSTRACT 

An automatic elevator controller for the final phase of an 
Instrument Landing System approach is designed using optimization 
theory and the practicality of the controller investigated. The 
problem is discussed and the assumptions stated. Then a mathe- 
matical model for the aircraft and a desired flare-out approach 
path are derived. The aircraft and approach limitations are estab- 
lished and the model is testel®. DynJiic Programming and the 
Parametric Expansion Method provide the optimal control from which 
the controller is designed. A computer program is developed to 
investigate the controller. The results are inconclusive and a 
recommendation for further study is made. 
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TABLE OF SYMBOLS AND ABBREVIATIONS 



T (Superscript) • 


3 matrix transpose operator 


. (Superscript) « 


« first time derivative 


, , (Superscript) ' 


» second time derivative 


. . . (Superscript) > 


» third time derivative 


IV (Superscript) • 


■ fourth time derivative 


d (Subscript) « 


> desired value of the variable 


e (Subscript) ‘ 


* denotes an equilibrium condition 


s (Subscript) < 


> denotes an aircraft constant 


A ‘ 


• 

» matrix relating z(t) and z(t) 




> an element of A matrix 


■ 


■ constants of the desired approach 


a • 


B an arbitrary instant of time 


B » 


« matrix relating z(t) and u(t) 


‘■i ■ 


> an element of B matrix 


D « 


* an arbitrary distance from the end of the runway 


E = 


> minimum value of the error index, J. called the 
mlnlmiun error function 


F = 


* matrix of weighting factors for system errors 
at the final time 


■ 


* an element of F matrix 


£(a) . 


■ an arbitrary function of time evaluated at the 
instant a 


h(t) * 


* aircraft altitude 


ILS » 


° Instrument Landing System 


J * 


» system error index 




* aircraft short period gain 
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k (t) = 

tap 


! time functions of E 


m s 


summation index 


P * 


summation index 


P. . = 

ij 


time variant coefficients of the system 
signal x^(t)Xj(t) 


Q(t) = 


time variant matrix of weighting factors for 
the system errors 


"ij “ 


an element of Q(t) matrix 


s = 


complex variable 


s 


aircraft path time constant 


t s 


' real time 


h = 


time aircraft is over end of runway on the 
desired approach 




' desired touchdown time 


u(t) = 


' elevator control signal 


V = 


aircraft velocity 


w = 


aircraft weight 


w 

s 


' aircraft short period resonant frequency 


X = 


aircraft x axis 


x(t) = 


measured signal vector 




' element of x(t) vector 


Z = 


aircraft z axis 


z(t) = 


' vector of changes in the measured signals 


^i(t) = 


element of _z(t) vector 


o^(t) = 


aircraft angle of attack 


“s “ 


angle of attack when aircraft stalls 



A B denotes changes in the variable about some 
equilibrium condition 
6(t) “ elevator deflection 

y(t) * glide path angle 
e(t) ■ pitch angle 

^ » aircraft short period damping factor 
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1. INTRODUCTION 



This study Is an attempt to sjmtheslze an automatic control 
mechanism for the terminal phase of an aircraft landing using 
optimization theory. The problem was originally presented by 
F. J. Ellert and C. W. Merrlan [2]. Dynamic Programming [1] and 
the Parametric Expansion Method [2] are used to design a simple 
controller. 

Section two describes the problem and presents the assumptions 
made to simplify the initial study. Originally the investigation of 
Ellert and Merriam was to be continued but after several weeks of 
investigation, it was found that their equations for the aircraft 
represented only the changes about an equilibrium condition [3] 
rather than the true signals. Equations representing only the changes 
are derived in Appendix A and are used in section three to obtain 
different equations for the actual aircraft signals. 

In section four a more realistic desired approach is defined 
and equations for the desired values of the measured signals are 
derived. The performance limits and system constraints are dis- 
cussed in section five. Section six covers the testing of the 
aircraft model. 

Section seven develops the theory to obtain the optimal 
control in terms of the measured signals, aircraft constants, and func- 
tions of time* The. necessary fiyst-order, linear differential equation 
to obtain the time functions for the optimal control are derived 
in Appendix B. Section eight outlines the mechanization of a 
controller using the optimal control equation. Section nine 
discusses the computer programming and the investigation 
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conducted with the results in section ten and the conclusions 



in section eleven. 
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2. PROBLEM STATEMENT 



Landing an aircraft In Inclement weather when the pilot has 
no visual contact with the runway until the last few hundred feet 
Is a function performed at many airports by an Instrument Landing 
System, commonly referred to as ILS« 

The aircraft follows an electronic beam which Is set to a 
glide path of approximately three degrees and provides an azimuth 
reference to the runway. As the aircraft descends, the beam resol- 
ution becomes poor and sizable altitude errors can occur. The three 
degree glide path Is normally followed to a point within half a mile 
from the approach end of the runway. The pilot then takes over 
visually to land the aircraft or Is waved off to try again. 

At this time, a choice must be made by the pilot to follow 
ILS and hold the rate of descent he has or to make a flare-out 
approach to touch down on the runway with a lesser rate of descent. 
The flare-out approach puts less stress on the aircraft and Is a 
more comfortable landing. 

A controller Is needed to accomplish the flare-out approach so 
that If visibility Is restricted, this kind of landing may still be 
made. The following assumptions are made to simplify these first 
studies : 

a. The aircraft Is laterally aligned with the runway. Only 
errors In the vertical plane will be considered. 

b. At one-half mile from the end of the runway, the aircraft 
will be waved off If altitude, rate of descent, pitch angle, and 
pitch rate exceed specified limits. 

c. Wind effects will not be considered. 
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3. SYSTEM MODEL DERIVATIOH 



The first step in the synthesis of an aircraft control system 
is the development of a mathematical model relating the control 
(in this case, elevator deflection) to some measurable response 
variables. It is essential that any model chosen be tested to 
Insure that it follows all known responses of the real system. This 
will be discussed further in section six. 

Figure 1 defines the aircraft coordinates, angles, and positive 
direction of elevator deflection. Equations (3.1) and (3.2) have 
been derived to describe the motion of an aircraft [2] under the 
following restrictive conditions: 

a. The equations are linearized with the assumption that 
deviations from an equilibrium flight path are small. Only the 
changes about this equilibrium condition are represented by 
equations (3.1) and (3.2). 

b. The glide path angle, y , is small enough that the small 
angle approximations, siny®Y and cos Y = 1 can be used. This 
approximation is valid for the landing paths to be investigated. 

c. It is assumed that the aircraft velocity, V, is held 
constant e 

Note that should a controller be found that is feasible from 
this basic synthesisj> a more realistic model of the aircraft should 



be used for further studies 



Lift 




Aircraft 

Reference 

Horizontal 

Flight Path 
— —X 



Z W 



Figure 1. Definition of Aircraft Coordinates, Angles, and Elevator 
Deflection. 



(3.1) A0(s) 



(I.* + 1) 

+ 2{M « + W ^ 

^ s s 



A6 (s) 



(3.2) A(s) 



(TgS + 1) 



AS 



(s) 



• • 

where A6(s) « Laplace transform of the change in pitch rate, a 6 (t). 

A6 (s) a Laplace transform of the change in elevator 
deflection, A6 (t). 

Ahi[s) » Laplace transform of the change in vertical 
acceleration. Ah (t), 

K ^ ^ constants with different values dependent on 

s s s 

the type of aircraft. 

Equations (3.1) and (3.2) can be combined to provide a transfer 
function relating the change in vertical acceleration and elevator 
deflection. 
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(3.3) Ah(s) = 



K VW 
s s 



2 2 
3 + 2rW s 4- W 

^ £ S 



A6(s) 



which can be written as a fourth order linear differential equation 



IV •« • 

(3.4) Ah(t) + 2^W ^*(t) + W ^Ah(t) = K VW A6(t) 
s s s s 

relating the change in the control variable to changes in the 

response variables: change in altitude, A h(t); change in rate of 

• • • 

ascent, Ah(t); change in vertical acceleration, Ah(t); and change 

••• 

in jerk, Ah(t). 

All of these signals are not directly measurable. Signals 
that are readily measurable in the aircraft are desired. These are: 

a. Altitude, h(t), which can be measured from a barometric 
altimeter or radar altimeter, 

b. Rate of ascent, h(t), measured with a barometric rate meter. 

c. Pitch angle, 0(t). 

d. Pitch rate, 0(t), both measurable from gyros. 

These are expressed in (3.5) as an equilibrium condition plus a 
change about that signal. 



(3.5) *,(t) = h (t) + Ah(t) 

X2(t) ® + Ah(t) 

Xs(t) = e^(t) + A0 (t) 

x^(t) = 0g(£) + A0 (t) 

where the subscript e denotes the equilibrium condition. 

A set of first-order differential equations can be developed 
from (3,2) and (3.4) for the changes about equilibrium and result 
in equations of the form 

(3.6) z{t) “ A^(t) + Bu(t) 
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where z(t) is a 4x1 vector and 

z^(t) = [Ah(t), Ah(t), A6(t), A0(t)] 

A is a known 4x4 matrix of aircraft constants 

B is a known 4x1 distribution vector of the scalar control, 

u(t). 



The derivation of (3.6) is given in appendix A and represents 
the change of the measurable signals, x(t), about the equilibrium 
condition. Equations (3.5) can now be written 
(3.7) Xj^(t) « hg(t) + z^(t) 



* 2 (t) = h^(t) + z^(t) 

«3(t) “ a^(t) + Z2(t) 

« 4 (t) - 6g(t) + z^(t) 

The equlllbrlvim condition, Xg(h)» is chosen to be a constant rate 
of descent, as variations about an ideal ILS approach are to be invest' 
Igated. x^(t) becomes 
(3.8) h^(t) = h(0) + fi(0)t 



hg(t) = h(0) 

6g(t) - 0(0) 

6g(t) = 0 

Then (3.7) can be expressed as 



(3.9) x^(t) - h(0) 


+ 


h(0)t + Z2(t) » Xj^(O) + X2(0)t + Zj(t) 


X2(t) - h(0) 


+ 


Z2(t) = X2(0) + Z2(t) 


X3(t) = 0(0) 


+ 


Z3(t) = x^(0) + Z3(t) 




where z^^(O) = Z 2 ( 0 ) = 


^3' 


(0) = 0 and z^(0) = x^(0). 
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A set of linear differential equations describing x(t) is 



(3.10) 


Xj(t) - 


6(0) + Zj^(t) 






X2(t) = 








i3<t) - 


^(t) 






=* 






substituting (A. 9; 


1 in (3.10) gives 




(3,11) 


X^(t) = 


h(0) + ^ 2 ^^^ 






*2(t) =* 








X3(t) = 


^4(^) 






^4(t) = 




+ a 44 Z 4 (t) + b^u(t) 


using 


(3.9), (3.11) can be written 




(3.12) 


Xj^(t) » 


X2(t) 






X2(t) “ 


322^*2^*^^ - X2(0)] + 


®23^^3^^^ - X3(0)] 




X3(t) = 


X4(t) 






x^(t) = 


a42[*2(t) - ^ 2 ( 0 )] + 


a43[X3(t) - X3(0)] 



which are the desired linear equations as a function of the signals 
x(t), and the control, u(t). The system model is also dependent on 
the equilibrium values of the descent rate and pitch angle because 
of the chosen equilibrium condition. 






4. DESIRED FLARE -OUT APPROACH 



At time t s 0 with the aircraft some horizontal distance, D, 
from the end of the runway, it is desired to have the aircraft leave 
the three degree glldeslope, increase its rate of descent, then flare- 
out to arrive over the end of the runway at time, t = t^, at a lower 
altitude than if it had remained on the glldeslope beam and with a 
lower rate of descent. From time, t = t^, until time of touchdown, 
t - t^, the rate of descent is held constant and all accelerations 
are zero. Figure 2 shows this desired path. An exponential altitude 
path over the time interval, < t < t^, would be more realistic 
but the constant rate of descent was chosen for simplicity. 




TIME (sec) 

Figure 2. Desired Flare -out Approach 
Note: All desired approach variables will be subscripted, d, to 

differentiate from the actual aircraft model signals. 

The following assumptions are made: 

a. Origin of the ILS glldeslope beam and the desired touchdown 
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point coincide 1000 feet down the runway. 



b. Prior to t = 0, the aircraft has been following an ideal 
ILS approach. 

From the problem definition and the general desired approach 
specif IcationSj, the known conditions must be compiled to derive 
equations for the desired measurable signals. 

First the known conditions at times t = 0, t « t^^, and t « t^ 
will be stated. Using these boundary conditions, equations will 
then be derived for the time intervals 0 < t < tj^ and t^^ < t < t^. 

Prior to time t => 0, the aircraft is descending on an ideal ILS 
approach. Since all measurable signals are continuous, the following 
conditions exist at t « 0: 

h^(0), h^(0), 0^(0), a^(0), Yjj(O) are constants and 
h^(0)> zero. 

The flight path specified over the interval t^^ < t < t^ is 
the same as before t = 0 (constant rate of descent). Continuity 
requires the same conditions at t^ and t^ with different constants. 

The desired flare-out approach path over 0 < t < in terms 
of altitude, h^(t), and derivatives can be approximated by a 
seventh order polynomial since there are eight boundary conditions, 
four at t =* 0 and four at t = tj^. These are 




• 3 4 5 6 

=■ + 5a^t + 6a^t + 7a^t 

2 3 4 5 

= 12a, t + ZOa^t' + 30a, t + 42a_t 







2 3 4 

= 24a, t + 60a-t + 120a,t + 210a_t 

4 5 6 7 
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where » [20 y( 0) + 15Y(t^)] V - 35 rhj(O) - h^(t^)] 

(t^)^ 

aj = [45 y( 0) + 39Y(tj)] - V + 84 th^j(O) - h^(t^)] 

ag = [36 y( 0) + 34Y(tj)] V - 70 [h^(0) - h^(tp] 

- [10y(0) + lOY(tj^)] ■ V + 20 fh^(O) - h^(t^)] 

(tp^ 

For the Interval < t < t^, the desired flight path is a 
constant rate of descent. The desired equations are: 

( 4 . 2 ) h^(t) » h^(tp + h^(tp[t-tj] 

hd(t) = h^(tp 

hd<t) » 0 

fr^(t) » 0 

The desired pitch angle, 6^(t), can be expressed as its equili- 
brium value, 0^(0), plus changes about 6^(0): 



(4.3) e^(t) - 6^(0) + A0 ^(t). 

Taking the derivative gives 

(4.4) 0^(t) = A0^(t). 

From (3.2), 

(4.5) A0 (t) = _Ts_ A”(t) + 1 A^r (t). 

° V V 

Integrating (4.5) 

(4.6) 40 (t) » Ts A h*(t) + 1 A h,(t). 

0 V V 

The desired values of (3.5) can now be written in terms of the 
derived equations 
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( 4 . 7 ) 



h^(0) + h^(0)t + 



x^^(t) 



0 < t < t, 



h^(tp + h^(t^)[t-«jl,,tj^ < t < t^ 

hd(0) + 4a^t^ + Sa^t^ + 6a^t^ + 7ayt^, 



X2d(t) 



< 



0 < t < t. 



< t < t^ 



^3j(^) 



■*■ (12a^t^ + 20a^t^ + 30a^t^ + 42a^t^) 

/ c ti 

^ + 1 (4a^t + Sa^t + 6a^t + 7a^t ) , 0 < t < tj. 



0jj(tj^), < t < tj 



^ (24a^t + 60a^t^ + 120a^t^ + 210a^t^) 

+ 1 (12a, + 20a_t^ + 30a,. t^ + 42a_t^), 

~ ^ 5 6 1 '* 

0 < t < 

0, t- < t < t, 

1 - “ f 



^4d(*'> " A 
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5. PERFORMANCE LIMITS AND SYSTEM CONSTRAINTS 



Aircraft constants, numerical values of the desired flight 
path, and variable limits must be specified before testing the 
system model. The T-28 aircraft was chosen for this Investigation. 
Values of the aircraft constants In (3.1) and (3.2) for the T-28 
In a landing configuration are: 



(5.1) = 


-1.350 sec”^ 


T 

s 


2.205 sec 


W = 
s 


0.6161 radlans/sec 


C = 


0.4230 


V 


157.8 ft. /sec 



The constants for equations (A. 7) are calculated from (5.1) 
with reference to figure 2. ILS weather mlnimums of visibility 
vary from a half to a quarter mile. Using this as a basis, a value 
of one-half mile Is chosen for D. In consistent units, 

(5.2) D = 2640 ft. 

Then from figure 2, 

(5.3) h^(0) = (D + 1000)tan 3° =* 190.8 ft. 

Using the small angle assumption, 

(5.4) h^(0) - y^(0)V 



where Y^(0) = 


-0.05236 radians 



and 



(5.5) t^ = 


D sec 
V 




(D + 1000) sec 
V 
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A value of twenty feet is chosen for the desired altitude 
over the end of the runway, then 

(5.6) ° 

® “0.02 radians 

d 1 

hd(ti> a Y^(tj^)V a “3.156 ft. /sec 

The physical limitations of the aircraft must be considered; 
important details are discussed below. 

1) Angle of Attack. The T-28 aircraft stalls at an angle of 
attack, a(t), of approximately 21 ° and a speed of 72 knots. As 

the angle of attack approaches this limit, the speed cannot be maln“ 
talned constant and the linearized system equations no longer repre- 
sent the aircraft motion. It is assumed that realistic representation 

o 

Is lost for a(t) > 18 . This value will be called the stall angle of 

attack, a . For the T-28 In landing configuration descending on an 

s 

ILS glideslope at 93.5 knots, the angle of attack, a(0), is 11° or 
0.1920 radians. 

2) Pitch Angle. The pitch angle at touchdown must be between 
the limits 0° < 0(t^) < 14°. The lower limit prevents the nose 
wheel from touching first and the upper limit prevents the tall skag 
from dragging on the runway. 

3) Elevator Control. The elevator motion is limited by the 
mechanical stops of the actuator arm. Assumed values are -35° and 
+15°. The direction is defined in figure 1. Aircraft motion with 
the elevator against these stops for non-zero time is not permitted 
for the linearized model used. 

From figure 1, 

(5.7) 0(t) ® a(t) + Y(t) 
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and the desired pitch angle at t » 0 Is 

(5.8) ^ ** radians. 

Since the angle of attack must remain below the stall value, 

01 , an equation for a(t) Is desired. From (3.2), the change In 
s 

pitch rate Is 

(5.9) A0(t) « Ts Ah(t) + __1_ ii(t). 

V V 

(5.7) may be written In the form 

(5.10) 6 (t) + A9(t) - a^(t) + Aa(t) + y (t) + Ay(t). 

e e c 

Taking the time derivative with 9g(t) = ag(t) = Yg(t) “ 0, 

(5.11) A*e(t) » Aa(t) + AY(t) 

and equating (5.9) and (5.11) gives 

(5.12) Aa(t) + Ay(t) = Ts Ah(t) + JL. Ah(t). 

V V 

It can be shown [3] that Aa(t) » T Ay so the T term of (5.12) Is 

s s 

associated with Aa(t) and 

(5.13) Aa(t) - ^s Ah(t) 

V 

(5.14) Ay(t) - _l_ Ah(t). 

V 

Integrating equations (5.13) and (5.14) yields 

(5.15) a(t) o a (t) + _£s Ah(t) 

® V 

(5.16) y(t) “ Y.(t) + 1 Ah(t) 

® V 

where a^(t) ■ a(0) and y^(t) = y(0). 

From (5.4), y (0) *» h(0) and (5.16) becomes 

V 

(5.17) y(t) = x,(t) 

V 

Using (5.7), ot(t) can now be determined from the system model 
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by 



(5.18) a(t) » XjCt) - _1_ x^Ct) 

In section 2, It was stated that 1£ the measurable signals » 
x(o), were outside specified limits, the aircraft would be waved 
off for another approach. These limits are somewhat arbitrary and 
should be chosen only after a close study of existing conditions of 
actual aircraft making ILS approaches. 

For this study the following limits are assumed: 



(5.19) x^(0) 




+ 


20 ft. 


^2^0) 






2 ft. /sec 


X3(0) 


■ ’‘3d<“) 


+ 


0.0350 radians 


*4<0) 


- *4d<“> 


+ 


0.0043 radians 
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6. MODEL TESTING 



Every mathematical model derived for a physical system should 
be tested by applying known controls to Insure that the model 
accurately describes the system. Originally the ^(t) equations were 
used as the aircraft model assuming that they represented the 
signals, problem there are two controls that can be 

used for testing. The control, u(t) « 0, applied to the model allows 
the aircraft to follow its equilibrium flight path which Is a constant 
rate of descent to the runway. From (3.4) a control can be derived 
using the equations, Xj(t), since they are continuous and known 



functions of time. 



(6.1) u(t) 






+ 



2€ A h (t) 

K VW 
s s 



Since ^g(t) and h^(t) equal zero. 



+ 1 A il,(t) 

K V 

s 



(6.2) = h^(t) 

Ahj(t) - h*(t) 



Using (4.1) and (4.2), 



(6.3) i^(t) - 


12a, t^ + 
4 


20a^t^ + 


30a^t^ 


5 

+ 40a^t 




24a^t + 


2 

60a^t + 


120agt^ 


+ 210a^t^ 


Ah^(t) - 


24a^ + 


120a^t + 


360agt^ 


+ 840a^t^ 


for 0 < t < tj^. 

(6.4) Ahj(t) - 


Ah*^(t) - 


AhY(t) = 


' 0 





for < t < t^ , 

With (6.3) and (6.4) used in (6.1) over their respective time intervals 
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and this u(t) applied to the model, it should follow the desired 
flight path shown In figure 2. 

When these controls were applied to the first model. It did 
not follow either flight path (see figure 3), After deriving the 
model In section 3, the same two controls were applied and the 
proper flight paths were followed (see figure 4). 
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7. APPLICATION OF OPTIMIZATION THEORY 



The first step in applying Dynamic Programming is to obtain 
an error index. This index determines the form of the resulting 
optimal control. The errors, x(^) " 0 < a < t^; 

the errors at the final time, x(t^) ~ Xj(t^); and the control, 
u(t), are considered to be significant for this problem. The Bolza 
form of error index is considered adequate to Insure satisfactory 
system performance. The error index, J, is defined by 

(7.1) J =« 1 |[x(tf) - x^Ctp]"^ F[x(tp - Xj(tpjj 

2/ o *- 

where F is a 4x4, positive, seml-deflnlte, constant matrix. 

Q(t) is a 4x4, positive, semi-definite, time variant matrix. 

Both F and Q(t) are arbitrary weighting functions of the 
system errors and are chosen to be symmetric. 

The condition by which (7.1) is a minimum is given by [1] 

(7.2) min [j. J" [x(a) - x,(a)]^Q(a) [x(a) - x,(a)] + u^(a)l+ ^ j » 0 

u(a)|_2 I J daj 

where a is an arbitrary instant of time in the interval,© < a < t^. 

t^ is the fixed final time. E is defined as the minimum error 

function, E[f(a), x(a)J = min J, which is assumed to be only a 

u(a) 

function of a and the signals measured at a. 

As a — ► tj, It can be shown that 

(7.3) E[f(t^), x(tpj = i|(x(tp - x^(tp]’^F[x(tp - 

Define 

(7.4) H(x(a), u(a)] = l|^[x(a) - Xd(®)l^Q^®n 2 (a) - x^(a)] + u^(a)| 
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then (7.2) can be written as 

(7.5) min |^[x(a), u(a)] + dE l 
u(a) L daj 



= 0 



The derivative of E with respect to a may be written as 

(7.6) ^ = d Eff(a). x(a)1 

da da 

4 

= df(a) + T" dxn ^Eff(a), x(a)l 

^f(a) da n=l da ^x^(a) 

The minimal error condition (7.2) can then be written 

(7.7) min rH[x(a), u(a)J + bE df(a) + x^ (a)bE + 

u(a) U( ' ’ ‘ 



^i£(a) da 



ix- 



Xo(a) + Xo(a) + x^(a) = 0 

^X2 ^x^ J 



For (7.7) to be a minimum for all a, 0 < a < t^, the first 
derivative of the term in brackets with respect to u(a) must be zero 
and the second derivative must be positive. E was defined as a 
function of a and x(a) and is not a function of u(a). From the 
system equations (3.12) it is seen that x^(a) is the only time 
derivative that is a function of u(a). The derivatives of (7.7) 
with respect to u(a) can then be written 
(7.8) XH + hE 'hh(a) = 0 



'bu(a) 

bu(a)' 



+ 



From (3.12) and (7.4) 
(7.9) b^4(a) 



"buXa) 






^u(a) 

bu(a)^ 



b^^4(a) 

^>u(a)^ 



> 0 



M 



'6u(a) 



u(a), ^ H 

^ u(a)^ 



= 0 



= 1 
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Substitution of (7.9) in (7.8) yields 

(7.10) u(a) = -b^ 'h E 



1 > 0 



(7.10) defines the optimal control for this system over the interval, 



0 < a < t_. Substituting (7.4), (3.12), and (7.10) in (7.7) gives 

* ^ 2 

(7.11) i |[x(a) - x^(a))^^ Q(a)[x(a) - 



id 



2 

(a)U+ V "1 + 

^ K.J 



df(a) + X (a) + Va^^[xAa) - x-(0)] + 

of(a) da bx, L ^ 



^23^*3^®^ - * 3 ( 0 )]! AE.+ 3c^(a) + [^421^*2^®^ " 

+ a^3[X3(a) - X3(0)] + a^^x^(a)]^ - = 0 

L^J 



which is a partial differential equation defining E. If (7.11) can 
be solved for E, then the optimal control of (7.10) can be found. 

The Parametric Expansion Method [2] is now applied to solve 

(7.11) and provide an optimal solution. The method is outlined as 
follows: 

Assume a fom for E consistent with the functional description. 

4 4 4 

(7.12) E a k(a)" - 2 2 ^ (a) + S 2i k (a)x (a)x (a) 

m=l m®! p«*l ^ ^ 

where k (a) ® k (a), 

mp pm 

Taking the derivative with respect to x^(a), substitution in (7.10) 
yields the optimal control as a function of time, the aircraft para- 
meters, and the measurable signals, x(^)« 

(7.13) u(t) a 2b^[k^(t) - k^^(t)x^(t) - k 2 ^(t)x 2 (t) 
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since a Is an arbitrary Instant, It can be real time, t. Obtain 
the derivatives required by (7.11) from (7.12). Substitute these 
In (7.11) and collect terms so that the expression Is In the form: 

(7.14) + P^Xj^(a) + P2X2(a) + P2X^(a) + P^x^(a) + Pj^j^x^^(a) 

Pl3Xl(a)x3(a) + Pj^^Xj^(a)x^(a) + P23X2(a)Xj(a) + 

^ 24*2 + ^’34’'3<®)x4(®) " 0 

The P coefficients contain terms of the k-parameters, time 
derivatives of the k(a)'s, the desired signals, the weighting factors, 
the aircraft constants and the chosen equlllbrltim conditions. For 

(7.14) to be valid for all values of the measured signals, each P 
coefficient must Independently equal zero. 

(7.15) = 0 

o m mp 

where m = 1, 2, 3, 4; p =» 1, 2, 3, 4. 

Equation (7.15) gives a set of 15 Independent, ordinary, first- 
order differential equations which define the k(a)*s. These equations 
with the boundary conditions are derived In Appendix B. 
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8. IMPLEMENTATION OF A CONTROLLER 



Several types of controllers could be designed depending on 
the asstuned form of E and the methods available for measuring time 
and velocity. The block diagram shown In figure 5 Is obtained from 
(7.13) and Is a simple analog control system. A constant speed motor 
rotates the shaft at a constant velocity. The k(t)'s of (7.13) 
represent the settings of wire-wound potentiometers where the wiper 
arms rotate with the shaft. The measured signals are multiplied 
by the appropriate k(t), summed and multiplied to provide the ele- 
vator deflection signal. 




Figure 5. Elevator Controller 

In the derivation of the k(t) equations, the equilibrium values 
of the descent rate, X 2 g(t), and pitch angle, ^ must be known. 
This means that the optimal k(t) equations are dependent on the 
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equilibrium conditions. This Is a function of the aircraft model 
and would be altered with a different aircraft model. Little error 
should be Introduced If these values are assumed to be the desired 
descent rate and pitch angle at problem start. 

(8.1) = ^2d^°^ 

By assuming these values the system Is no longer optimal, but It Is 
desired to sacrifice optimality In the Interest of simplicity and to 
Investigate the feasibility of this simple controller. The k(t) 
equations are also dependent on the aircraft constants so a diff- 
erent set must be obtained for each type of aircraft. 

At this point two big questions remain. 1) Can k^(t), kj^^(t), 
k 2 ^(t), k^^Ct), and k^^(t) be found such that, for all aircraft 
with x(0) within the signal limits given by (5.19), the aircraft 
will make a satisfactory landing? 2) Are the k(t)'s found above 
linear enough to be wound on potentiometers? A method to investi- 
gate the answers to these questions Is outlined In section 9. 
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9. INVESTIGATION PROCEDURE 



Since all of the constants and variables of the fe(t) equations, 
(B.3), have been calculated or are known except the weighting factors 
F and Q(t), the answers to the questions in section 8 can be investi- 
gated by the following procedure; 

1) Arbitrarily pick the weighting factors of the error index, 

F and Q(t). 

2) Integrate equations (B.5) from t = t^ to t = 0 to obtain 
the k(t)'s of (7,13). 

3) Integrate the aircraft model equations (3.12) from 

t ■ 0 to t " tj using u(t) of (7.13) with the k(t)'s obtained in 2, 

4) Observe the final signal errors, x(t.) - x,(t_). If these 

~ r ""d t 

are not within prescribed limits, change F and/or Q(t) and repeat 
steps 2 through 4. 

A program was written to accomplish the above procedure using the 
CDC 1604 computer and an existing modified Runge-Kutta integration 
subroutine. There were two main sections in the computer program. 

1) Equations (B.5) were Integrated from t = t^^ to t = 0 with the 
values of k^(t), kl4(t)» »^24<‘^>' k 34 <t), and k 44 (t) stored in a 
buffer. 2) Three sets of equations (3,12) were Integrated from 

t =« 0 to t “ tj using u(t) of (7.13) and the k(t)'s that were 

stored in the buffer. 

Three aircraft models were chosen for each program so that one 
could be started with positive errors, one started on the desired 
trajectory, and the third started with negative errors. This allowed 
all of the system errors to be observed in one program. 
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The final error limits were arbitrarily chosen as: 

a) An altitude, descent rate combination such that all three 
aircraft land within +500 feet of the desired touchdown point with 
a descent rate of + 100 ft. /min or + 1.66 ft. /sec from the desired 
rate of descent, *2d^^P* 

b) Pitch angle as stated in section 5, 0° < 6(tj.) < 1A°. 

c) No limit was set on the pitch rate but it must be small 
for the approach to closely follow the desired flight path. 

For this investigation, the aircraft was considered to have 
only an altitude error and Q(t) = Q, a constant matrix consisting 
of the diagonal elements only. The study was conducted in the 
following manner: 

1) Determine if successful approaches can be made with con- 
stant weighting of any one signal; or any two signals; or combin- 
ations of any three signals; and finally, combinations of all 
four signals. 

2) Determine if successful approaches can be made with constant 
weighting of all signals, progressively adding diagonal elements of 
the F matrix. 

3) When successful approaches are found, determine the nature 
of the k(t)'s required for the controller. 
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10. RESULTS 



Digital simulation was used to study the areas outlined in 
section 9. The results were: 

1) Constant weighting of all four signals was required to 

obtain a satisfactory trajectory. All final error limits were 
met except that of the minimum descent rate. Figure 6 shows the 
altitude vs. time/distance for three aircraft started with an 
initial altitude error of + 10 feet, where = 2.5xl0’^^, 

q22 “ 8.5x10 q^^ “ 6.0x10*^, and q^^ =« 8.8xl0~^. All 

other q's and F were zero. The minimum descent rate at touchdown 
was -5.36 ft. /sec compared with the limit of -4.91 ft. /sec. 

2) This area was not fully investigated but successful 

approaches were made with constant weighting of the four signals 
(Q equal to zero except q22» 933* 9^^) weighting all 

four signals at the final time (F equal to zero except ^22* 

fj^j ^44)* Figure 7 shows the altitude vs. time/distance for three 
aircraft started with an initial altitude error of + 20 feet where 



‘111 = 


7.10x10"^? q22 » 


2.95x10“^, q^j = 7.00x10"^, 




-4 


-9 « -7 


944 = 


8.00x10 and fj^j^ 


= 1.80x10 f22 = 2.00x10 , 


%3 “ 


2.90x10"^, = 


.3 

3.00x10 . All final errors were within 



specified limits. 

3) For the approach of part 1, ^rtiere F is zero, k^(t) is shown 
in figure 8, k^^Ct) in figure 9, and k^(t) in figure 10. 
k2^(t) were insignificant- (i.e. nearly zero). 

For the approach of part 2, k^(t) is shown in figure 11, 
in figure 12, and k^(t) in figure 13. Again kj^^(t) and k^^Ct) were 
insignificant . 
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11. CONCLUSIONS 



For the results specified in section 10, the following conclu- 
sions are stated: 

1) Although the minimum descent rate error was not within 
limits using the Q matrix only, it is felt that with several small 
adjustments of the weighting factors the descent rate could be brought 
within limits. The k(t)'s are of the same sign with maximum values 
less than one and could be wound on the potentiometers. 

2) While all of the error limits were met for the approach 
using both F and Q, the k(t)'s changed sign and are unsatisfactory 
for this controller. As adjustments were made to F and Q, bringing 
the flight paths closer to the desired trajectory, the k(t)'s 
approached the same sign over the Interval of flight. It is believed 
that if F and Q are adjusted for flight paths with smaller errors, 
this problem would be eliminated. 

Two recommendations are made: 

1) That investigation be continued until definite conclusions 
can be obtained about the practicality of this controller. 

2) That a method be investigated to systemmatically adjust 
the weighting factors, the elements of the Q and F matrices. 
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APPENDIX A. DERIVATION OF £(t) EQUATIONS 
Define 



(A.l) z^(t) « Ah(£) 






Z3<t) 

z^it) 



th(t) 

A9(t) 

A0(t) 
u(t) = Lg(t) 

The time derivatives of the z(t)'s are: 
(A. 2) Zj^(t) « Ah(t) = z^it) 

d£(t) 

A0(t) 

A*9(t) 



Z2<t) 

Z3<t) 



Z4<t) 



z^Ct) 



Equations for Ah(t) and A9(t) in terms of the states, £(t), are 
needed. (3.2) can be written as 

VA9(s) 



(A. 3) T Ah(s) + A^(s) 

s 



Taking the inverse Laplace transform and solving for Ah(t) yields 

(A. 4) Ah(t) = -_iAh(t) + _V^0(t) 

T T 



Differentiating (A. 4) gives 

15^ • 

(A. 5) Ah (t) = -JLAh(t) + _V<A9(t) 

T T 

s s 

• • • 

Integrating (A. 4) and assuming Ah(0) = 0 gives 

(A. 6) 2ih(t) = -_lAh(t) + _VA9(t) 

T T 

s s 

Substituting equations (A. 4), (A. 5) and (A. 6) in (3.2) and solving 
for A9(t) gives 
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(A. 7 ) A 0 (t) = 1 

V 






1 - 2SWs 

™ 2 T 
T s 

s 



- 2S'Ws 



+ w ae(t) + r_i 

J 



2^- 



- 2SwJ 4e(t) 



+ Ws^s<‘> 



Inserting (A. 6) and (A. 7 ), written in terns of (A.l), into (A. 2 ) 
yields 

(A. 8) 2^(t) =■ z^(t) 

Zo(t) “ - _1 z (t) + _JV z,(t) 



2j(t) = z^(t) 



^ 4 (t) 



22^*=) • 



which is in the form 

(A. 9 ) z(t) “ Az(t) +Bu(t) 

where 



A s 



1 fjL - 2 |w£ + w 
V [t/ ^s 

+ - 2?W^j z^(t) + KT 



1 - 2SWfi + W 



Z3it; 



'^s 



8^8 









0 10 0 




'0 


0 a^2 S 23 0 


B s 


0 




1 




0 0 0 1 




0 


0 ^42 ®43 ®44 







and a22 » 323 - “ | 

s 



‘43 



2^8 + W 



T T 

s s 



1 - 2 t^a + W ^ 

— 2 ® 


• ®44 “ 


1 - 2lW 
f" ® 


T 

_ s s J 







, b,»K T W 
* 4 s s s 
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let 



APPENDIX B. DERIVATION OF k (t) EQUATIONS 

mp 

The partial derivatives of E defined in (7.12) are needed for 
substitution in (7.11). Since a is an arbitrary Instant of time in 
the interval 0 < a < t^, it can be chosen as real time, t. The 
functional notation, (t), on all the variables will be dropped for 
simplicity. 

The required partial derivatives of E are: 

(B.l) ^E df(t) = )jE k + ^ kj^ + ^E L 



'bf(t) dt 






ik. 






+ ^ ^ \ ^11 + kj3 



ik. 






11 



22 



^k 



33 



+ ^E k^^ + >jE ^ E k^^ + b E k 



i)k 



44 






12 



12 



^k 



13 



ik 



14 



14 



ik 



23 



23 



^ ^24 ■*" ^ ^4, 



24 



ik 



34 



4e = 1,11= "2Xi>lE = -2x , iE = -2x., 

Tk >k. ^ ik- ^ >k 



ll = 
\k^ 


“2x^ , 


'bE a 

Uu 


^22 


2 

Jk33 


'bE » 

^k44 


2 

^4 


'^E a 
^2 


2X,X , ^ a 

^3 


2x X , bE_ = 2x X. , 


a 


2x«x 


■a, ^E 


a 2x„x, , ^E = 


2x,x. , 


^3 


2 


®‘^24 


^ ^ w 

4^34 


3 4 



bE 

lE_ 



^1 ^ ^11*1 ^ ^12^2 ^13*3 ^14*4^* 



2(*k2 "f* ^12^1 ^ ^22^2 ^ ^23^3 ^24^4^* 



2(*k3 “I" ^2.3^1 ^23^2 ^ ^33^3 ^34^4^ * 



" 2(-k^ + + k2^X2 + + ^44*4)* 

4 




^ “ 4(k/ 



2 ^ 4 ^ 14*1 “ ^^ 4 ^^ 24^2 ”^^ 4 *^ 34^3 



- 2^44^ huW + ^24^^2^ + ^4^*3^+ ^^44^*4^ 
■*■ 2 ^ 14 ^ 243 ^ 1*2 ■*■ 2^14^34^1*3 ■*■ 2^14^44*1*4 

+ + 2k2^k^^^X2*^ + 2kj^k^^X2X^) 

The derivatives of (B.l) are then substituted in (7.11) to obtain 
(B.2) 1 j^Cx - x^)^Q(x - x^)| -2b^^(k^^ - 2k^k^^x^ “2k^lc2^X2 

”^'^4'^34*3 ■ ^^^^44*4 *^14 ^1 ^24 ^2 ■*■ *^34 ^3 

^^44 *4 '*' ^^14'^24^l’'2 '*' ^'^14^34*l’'3 '*’ 2'^i4'^44*i*4 
2^24^34*2*3 + 2k24\4'^2''4 + 2k3^k^^X3X^) + k - 2k^x^ - 
2k*2*2 - 2^3X3 - 2k^x^ + k^^Xj^ + k22*2^ + ^33*3^ + 

^4*4 ■*■ 2^12*1*2 2^i3Xj^X 3 + 2kj^^XjX^ + 2k23X2X3 + 

2^24*2*4 2^3^X3X^ + 2x2(-kj^ + + kj^2*2 ^13*3 

^ 14 * 4 ^ ■*■ 2 f ® 22^*2 “ * 2 ^°^^ ■*■ ® 23^*3 “ * 3(°^^1 ^"^2 ■** 

^12*1 '** ^^22*2 '*’ ^23*3 '*’ *^24*4^ 2x^(-k3 + k^^3Xj^ + 

^23*2 '*’ ^^33*3 '*' ^34*4^ '*’ ^^*42 ^^2 “ * 2 ^°^^ ''' 

®43^*3 “ *3^®)^ ^44*4! ^"^4 '*’ ^14*1 ^24*2 '*' 

^34*3 + ^4V “ ° 



The matrices in the first term are expanded and all terms are 
grouped in the form of (7,14). The resulting P coefficients are: 
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(B.4) P ® k - 2b/k/ + 
o 



4 4 



L*ld **“ *J22*2d ■'■ ‘*33*3d 



^44*4d ^ ■*■ ^I2^1d*2d *^13^1d*3d ■*■ ‘*14*ld\d *^23*2d*3d 

‘^ 24 "' 2 d'' 4 d + ‘» 34 '' 3 d^d + 2k2 (822*2(0) + 823X3(0)) + 2k^(a^2'^2(°> 



■*■ ®43^3 



( 0 )) 



?! e -2kj^ + 4b^ ^*^14 “ (^^ii^id ■’■ *^12^2d ‘^13*3d '*' 

<114*44) - 2>'l2(=22*2<®> >23'‘3<®>) - * 

* 43 * 3 ( 0 )) 

?2 » “2k2 + 4b^ ^4^24 " ^‘*12*ld '*' ^22^2d ‘^23*3d ■*■ 

<l24*4d) - 2«22(*22='2<“> + *23*3<®>) ‘ ^‘‘24(*42*2^“> + 



®43^3 



(0)) - 2k^ - 2822^2 - 28^2^4 



^3 ^ ”^*^3 ■'■ ^^4 \^34 ‘ ^‘•13^'ld ■*■ ‘^23*2d '*' *^33*3d ■*■ 

*l34*4d^ “ ^^23^22*2^°^ ■*■ ®23*3^°^ " ^^ 34 (^ 42 * 2 ^°^ 

®43*3^®^^ ° ^®23^2 “ ^®43^ 

^4 “ ^^4 V 44 " ^‘*14’'ld *^24*2d ^34*3d 

‘l44*4d^ *’ ^^24{®22*2^®^ **■ ®23*3^°^ “ ^*'44{®42*2^°^ ‘'’ 

( 0 )) 



®43*3 



2k^ ° 2a^^k^ 



^11 ^ ^11 '*' ”^11 “ ^14 



^22 ^ ^22 "'■ ■|%2 “ ^\^24 ‘‘' ^*^12 ■*■ ^® 22^22 ''' ^® 42’^2 



22 “'22 “ ■|'*22 ‘■"'4 “24 ’ ““12 ' ““22“22 ’ "''42*24 

P^-i ® k,» + - 2b, + 2a-,k-_ + 2a, „k, 



2, 2 



^44 ® *^44 ■*■ i%4 " ^44 **■ ^^34 ^®44^44 
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12 


a 


2^12 + ^12 ■ ^^4^14*^24 + 


2kii 4* ! 


13 


a 


^^13 *^13 " ^^4 ^14*^34 


^®23^12 


14 


a 


^^14 ■‘‘ *^14 ‘ ^^4 *^14^44 


2ki3 + : 


23 


a 


2^23 + q23 “ ^^4^*^24^34 


2ki3 + ! 






2®42*^34 ■*■ ^®43^24 




24 


- 


^^24 + ’24 - “>4S4'‘44 + 


2ki^ + : 






+ 2®42*^44'*‘ ^®44*^24 




34 


= 


2*^34 ■'■ *^34 " ^^4^^34^44 ■*■ 


2k33 + 






+ ^®44^34 




each P 


coefficient to equal zero 


for all 






measured signals as required by (7.15), the following first order, 
linear differential equations must be satisfied; 



(B.5) 

k = 2b^\^ - + q22*2d *^33*3d '^44*4d ^ 

‘*12*ld*2d " ‘*13^1d^3d “ ^14*ld’^4d ” '•23*2d*3d 
- ‘l24''2d''4d " ^34'^3d''4d “ 2^2(822*2(0) + «23*3(°>) 

* 2k4(a42X2(0) + a^3X3(0)j 

ki = 2b^\^kj^^ - + qi2*2d ^ ^13*3d ‘*14^4d^ 

“ ^12(®22*2(°^ '** ®23*3(°^) “ ^14(®42*2^°^ ®43^3(®^) 

^2 “ ^^4 *^4*^24 “ “^‘’l2*ld **22*2d ■*■ ‘*23*3d *^24*4d^ 

"^22(®22*2(®^ ®23*3(°^) " ^24C®42*2(°^ ®43’'3^°^ 

■ ‘‘l • ^22*^2 • ^2^ 
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S “ ^^4 ^4^34 “ ■*" ^23*2d '*' ‘*33^3d ^34*4d^ 

2 

“ *^23(^22*2^^^ ^ ®23^3^°^) * *^34C®42^'2^®^ ■*■ ®43*3^°^^ 

“ ®23*^2 “ ®43*^4 



^4 


a 


^^4^4*^44 * |^^14*; 


Ld “*■ ‘*24^2d **34^3d ‘*44*4d^ 
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where the boundary conditions at t ® t^ must be specified. 

Equation (7.3) is the expression for E at t « t^. Equating 

this to (7.12) evaluated at t = t^ gives 
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Expanding both sides and solving for the k(tp's gives the boundary 
conditions on the k's: 
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